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SUMMARY 

Stress distributions were calculated near a circular hole in [0/90]^ and 
[90/0]g laminates, using a three-dimensional finite-element analysis. These 
stress distributions were presented three ways: through the thickness at the 

hole boundary, along radial lines at the 0/90 and 90/0 interfaces, and around 
the hole at these interfaces. The interlaminar normal stress, a^, and the 
shear stress, o^q, distributions had very steep gradients near the hole 
boundary, suggesting interlaminar stress singularities. The largest compres- 
sive a, stress occurred at about 60° from the load axis; the largest a „ 

Z Zu 

occurred near 75°. A simple procedure was introduced to Calculate interlami- 
nar stresses near the hole boundary. It used stresses calculated by an exact 
two-dimensional analysis of a laminate with a hole as input to a quasi-three- 
dimensional model. It produced stress.es that agreed closely with those from 
the three-dimensional finite-element model. For laminates with holes, this 
simple procedure may become a viable alternative to three-dimensional finite- 
element analyses. 
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INTRODUCTION 

Composite laminates are being used extensively in aerospace structures. 
Many such laminates develop high interlaminar stresses near holes; such 
stresses can cause delamination. Unfortunately, the stress state near a cir- 
cular hole in a laminate is a complex three-dimensional (3-D) problem. The 
only exact solutions to the problems are two-dimensional (2-D) [1]. The 
literature contains little work on 3-D problems. Dana and Barker [2] used the 
finite-element method to analyze the stresses around circular holes in [0/90]^, 
[90/0]g, and [±45]^ laminates; they used very coarse models, however. Experi- 
ence with a less complex straight-edge problem [3] showed that coarse models 
may not yield accurate results near the interlaminar stress singularity that 
can exist at the hole boundary. Tang [4] used a boundary-layer theory to 
obtain stress distributions around a circular hole in [0/90]^ and [±45]^ lami- 
nates, but he used an extremely large value of 400 for the ratio of hole 
radius to ply thickness. 

One purpose of this paper is to present accurate 3-D stresses for a com- 
posite laminate with a more representative ratio of hole radius to the ply 
thickness. A second purpose is to introduce a method for estimating the major 
interlaminar stresses near a circular-hole. The present paper does not pre- 
sent design information; instead, it addresses the basic understanding of 3-D 
stresses at a laminate hole. 

The present study focused on four-ply laminates with [0/90]^ and [90/0]^ 
stacking sequences. A central circular hole of radius equal to 20 times the 
ply thickness was considered. The laminate had a length-to-width ratio of two, 
and the width was six times the hole diameter. A laminate of such dimensions 
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is, effectively, infinitely wide and long relative to the hole diameter. This 
laminate was analyzed by a finite-element method, like that in Ref. 3, with 
very high mesh refinement near the intersection of the interface between plies 
and the hole boundary. 

The stress distributions are presented in three ways; through the thick- 
ness, along radial lines away from the hole, and around the hole. Because an 
interlaminar stress singularity was expected between the 0° and 90° plies at 
the hole, the computed stresses were given very near but not at the hole 
boundary. Furthermore, the stress magnitudes very near this singularity may 
depend on finite-element mesh refinement. As a result, the peak stresses near 
the hole were used only for relative comparisons at different locations around 
the hole and for relative comparisons among the 3-D stress components at speci- 
fic points. 

SYMBOLS 

Young's modulus in i -direction, MPa 
G. . shear modulus, MPa 

* J 

h ply thickness, m 

L half-length of the laminate, m 

R radius of the circular hole, m 

r,e polar coordinates 

Sg gross-section stress, MPa 

u,v,w. displacements in x-, y-, and z-direction, respectively, m 

W half-width of the laminate, m 

x,y,z Cartesian coordinates, m 


Poisson’s ratios 
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{a} cylindrical stresses o^, o^, a^, MPa 

{e} cylindrical strains e^, Eq, e^g, Eq^, e^^ 

Subscripts 

i.j i.j = 1.2.3 

1,2,3 longitudinal, transverse, and thickness directions of a ply 

ANALYSIS 

This section describes the problem analyzed and the corresponding finite- 
element model . 


Laminated Composite Plate with a Circular Hole 

Fig. 1 shows a four-ply laminate with a half-length of 240h, half-width 
of 120h, and half-thickness of 2h. It has a central circular hole with a 
radius of 20h. The ends of the laminate were subjected to a uniform displace- 
ment in the x-direction. The gross-section stress necessary for this deforma- 
tion was Sg. 

As previously mentioned, two laminates, [0/90]^ and [90/0]^, were con- 
sidered in the analyses. Each ply was idealized as a homogeneous, elastic, 
orthotropic material with the same properties as in Refs. 3, 5, 6, and 7. 


= 138 GPa (20 x 10° psi) 

E22 = E33 = 14.5 GPa (2.1 X 10® psi) 


®12 “ ®23 " *"31 " (0-85 X 10° psi) 


^12 " ^13 " ^23 " 


The subscripts 1, 2, and 3 denote the longitudinal, transverse, and thick- 


ness directions, respectively, of an individual ply. 
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The analysis is aimed at obtaining the stress distribution in the ideal- 
ized laminate due to the above-mentioned loading, S^. Because of the com- 
plexity in satisfying the equations of equilibrium, compatibility conditions, 
and boundary conditions of orthotropic elasticity, exact solutions are not 
available. Therefore, approximate methods are being used to obtain the stress 
distributions in the laminate. In the present study, a finite-element analy- 
sis based on displacement formulation was used. 

Because of the symmetries in the problem, only one-eighth of the laminate 
was analyzed, as shown in Fig. 1. The displacements u, v, and w were pre- 
scribed as zero on x = 0, y = 0, and z = 0 planes, respectively. On the 
X = 240h plane, a constant displacement u = u^ was prescribed. The gross- 
section stress Sg necessary for this deformation was 302uQ/h, MPa. 

To facilititate presentation of results, a cylindrical coordinate system 
(r,6,z) was also used. The stresses in the cylindrical system are shown in 
Fig. 1. 

Finite-Element Idealization 

One-eighth of the laminate (0 < x < 240h, 0 < y < 120h, and 0 < z < 2h) 

was modeled by 20-node isoparametric brick elements. The plan view of the. 
model is shown in Fig. 2(a) and a 20-noded brick element is shown in Fig. 2(b). 

The model was subdivided into three regions: a 2-D region, a transition 

region, and a 3-D region. The unshaded region, in Fig. 2(a) corresponds to the 
2-D region. The shaded region very near the hole corresponds to the transi- 
tion and the 3-D regions. 

Fig. 3 shows the through-the-thickness modeling typical of the nine radial 
lines shown in Fig. 2(a). The 3-D region extended from 20h to 26h, the 
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transition was in the region 26h to 27h, and the 2-D region was from 27h 
onward. In the 3-D region very near the hole (20h < r < 21h and 0 < 2 < 2h), 
a polar-type arrangement of elements was used. The polar mesh had a very fine 
arrangement of elements very near the hole. This arrangement yields a concen- 
tration of stress values very near the intersection of the interface and the 
hole (z = h, r = R), where a stress singularity is expected. A very similar 
polar mesh was used in Ref. 3 for a quasi-three-dimensional analysis of the ‘ 
straight-edge problem in a laminate. 

The 2-D region extends outward from r = 27h (see Fig. 3). This 2-D 
region was also modeled by 20-node isoparametric brick elements. In this 
region, one element per ply was used through the thickness and the displace- 
ments were assumed to be independent of z. For example, nodes i, j, and k 
in Fig. 3 had identical u, v, and w displacements. Furthermore, the con- 
stitutive relations for the elements in the 2-D region were edited to produce 
plane stress (a^ = o) conditions. Thus, the 20-node isoparametric brick ele- 
ments degenerated to 8-node plane stress quadrilateral elements. 

The transition elements provided compatibility between the 3-D and 2-D 
regions. The transition elements were developed like those reported in 
Ref. 8. 

The element sizes used in the present finite-element modeling are an 
order of magnitude smaller than the fiber diameters used in current graphite/ 
epoxy laminates. Even so, the present analysis was not aimed at idealizing 
the fibers and matrix separately, but was aimed at analyzing laminates ideal- 
ized as consisting of homogeneous, orthotropic plies. Similar assumptions 
were also used in Refs. 2 and 4. 
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EVALUATION OF THE FINITE-ELEMENT MODEL 

The finite-element model for a 6 = constant plane is presented in 
Fig. 3. The polar arrangement of elements is much like that used in Ref. 3 
to analyze a straight-edge problem involving the same material properties and 
stacking sequences as the present problem. Refs. 3 and 7 show that such a 
polar mesh gives accurate stresses even near a singularity. Therefore, the 
present through-the-thickness modeling was expected to be accurate near the 
interlaminar singularity. 

The modeling around the hole was evaluated by comparing 2-D stress dis- 
tributions with a known exact solution [1]. Because of the polar arrangement 
of elements, the inplane modeling at the midplane (z = 0) and the top surface 
(z = 2h) was identical and was coarser than the inplane modeling at the inter- 
face (z = h) (see Fig. 3). The inplane modeling at the midplane was used in 
the 2-D finite-element analysis of a [0/90]^ ^or [90/0]^^ laminate with a 
hole. 

Fig. 4 shows the hoop stress distribution on the three radial 

lines, 0 = 0°, 45°, and 90°. The distance from the hole is nondimensional- 
ized with the ply thickness. The solid curves represent 2-D (midplane mesh) 
results and the dashed curves represent the exact 2-D solution for an infi- 
nite laminate [1]. The two solutions agree closely along all three radial 
lines. The exact hoop stress Og has a maximum at 6 = 90° and the corre- 
sponding stress-concentration factor, Kj, for an infinite laminate is 4.87. 
For the present finite-element model, the finite width and length elevate the 
Kj by about 2.4 percent [9]. Therefore, the exact Kj for the present 
model is about 4.99. The finite-element solution gave a Kj of 4.92, which 
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is only 1.4 percent below the exact value. Therefore, even the relatively 
coarse midplane modeling predicted the stress distributions accurately. 

The next parameter to be checked was the extent of the 3-D region. The 
quasi -3-D analyses in the literature [3, 5, 6, and 7] show that the 3-D 
stresses die away within a distance of about 4h from a straight free edge. 
For the present free-hole problem, two different 3-D models were considered. 
One model had 3-D elements extending to a distance 6h from the hole (from 
20h to 26h), as shown previously in Fig. 3. The second model had a smaller 
3-D region, extending to 4h (i.e., from 20h to 24h) from the hole. The 
3-D stresses from both models were nearly identical, showing that either model 
was adequate. The model with the larger 3-D region (20h < r < 26h) was used 
throughout the remainder of this study. This model had about 19,000 degrees 
of freedom. 

Because the finite-element idealization shown in Figs. 2 and 3 can 
account for the through-the-thickness variations, inplane variations, and 3-D 
effects, the model was considered accurate and was used to study the stress 
distribution in the [0/90] and [90/0] laminates with a central circular 

V ^ I 

hole. 

RESULTS AND DISCUSSION 

As previously mentioned, stresses are presented first as through-the- 
thickness distributions, next as radial distributions at the interface, and 
then as circumferential distributions at the interface. Lastly, estimated 
interlaminar stresses based on a procedure described in the Appendix are com- 
pared with stresses from the 3-D hole model. The results for the [0/90]^ and 
[90/0]g laminates are presented and discussed together. 
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Through-the-Thickness Distributions 

Hoop Stress, Og— Figs. 5(a) and 5(b) show the Og distributions through 

the thickness at the hole boundary (r/R = 1 ) for the [0/90] and [90/0] lami- 

s s 

nates, respectively. For comparison, stresses from the 2-D exact solution for 
an infinite plate [1] are shown as dashed lines. These 2-D stresses are, by 
assumption, uniform through each ply. The hoop stress, a., at 6 = 90° in 
the 0° plies of both laminates shows the largest variation. The gradient of 
Og in these plies is sharpest near the interface (z/h =1). This "peak" is 
probably caused by an interlaminar stress singularity at the interface [3]. 
Although the present finite-element method yields accurate stresses very near 
singularities, it cannot calculate stresses ^ such singularities [7]. For 
this reason, the stress results from the two elements closest to each side of 
the interface were neglected, as in Ref. 7. Because these elements were. within 
the narrow range 0.998 < z/h < 1.002, their absence is not noticeable in 
Fig. 5. (These two elements were also neglected in all subsequent through- the- 
thickness plots.) In general, the Og peaks for e = 0° and 45° in Fig. 5 
were quite small. Results for the two laminates were similar and also agreed 
reasonably well with the exact 2-D solution. 

Interlaminar Normal Stress, a^--Figs. 6(a) and 6(b) show through the 
thickness at the 9 = 0°, 45°, and 90° angular locations on the hole boundary 
at the [0/90]g laminate. At all three locations, the curves have steep gradi- 
ents near the interface. These results strongly suggest the existence of inter- 
laminar stress singularities at the hole boundary. Furthermore, at 9=0° and 
90° the stresses are discontinuous across the interface (z = h). Similar 
behavior was observed in Refs. 3 and 6. However, the present discontinuities 
may be somewhat magnified because the stresses were not plotted for the two 
elements on either side of the interface, as previously discussed. Near the 
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midplane the distributions have smooth gradients and show no evidence of 
a stress singularity. The interfacial is compressive at 0=0° and 45° 
and tensile at 0 = 90°. At the midplane, however, the is tensile at 
0 = 45° and compressive at 0=0° and 90°. 

Figs. 6(c) and 6(d) show o through the thickness for the [90/0] lami- 
nate.. As in the [0/90]^ case, the distributions near the interface show steep 
gradients and discontinuities, suggesting interlaminar stress singularities at 
the hole boundary. The interfacial is compressive at 0=0° and 45° and 
tensile at 0 = 90°, as in the [0/90]^ case, but with slightly different magni- 
tudes. At the midplane, in contrast to the [0/90]^ case, the is tensile 
only at 0 = 90°. 

Interlaminar Shear Stress, a^g— Figs. 7(a) and 7(b) present the a^g 
distributions through the thickness. Distributions are presented only for 
0 = 45° because a^g is identically zero at 0=0° and 90°, which are planes 
of symmetry. The two distributions are similar but of opposite sign. This can 
be explained because the load transfers from the 0° top ply to the midply in 
the [0/90]g case, while in the [90/0]^ case the load transfers from the 0° mid- 
ply to the top ply. Both laminates show steep gradients near the interface, 
again suggesting a stress singularity. 

Radial Distributions 

The radial distributions for interlaminar stresses a^, and 

are presented at three angular positions, 0 = 0°, 45°, and 90°. 

Interlaminar Normal Stress, o^— Figs. 8(a) and 8(b) show the interlaminar 
distributions away from the hole for [0/90]^ and [90/0]^ laminates, respec- 
tively. In these figures, and as in Fig; 4, the stresses are plotted against 
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(r - R)/h, the distance from the hole in terms of the ply thickness h. The 
interlaminar stresses from the two elements closest to the hole boundary were 
neglected because of the possible singularity there. Beyond these two ele- 
ments ((r - R)/h > 0.002), the discontinuities across the interface were 
small, only about 2 percent of that shown in Fig. 6(b) for the hole boundary. 
Therefore, the interlaminar values were calculated by averaging the 0° 
and 90° ply values at the interface, a customary finite-element practice. 

The distributions for the [0/90]^ laminate (Fig. 8(a)) have steep 
gradients very close to the hole boundary. The largest value of occurred 
at 0 = 45° and is about 15 percent of the gross section stress, S^. At all 
three locations, switches sign very near the hole boundary. Furthermore, 
for larger values of (r - R)/h, the becomes vanishingly small, with 
at 0=0° showing the most rapid decay to zero. 

Fig. 8(b) shows that the distributions for the [90/0]^ laminate are 
Somewhat similar to the [0/90]^ case. However, in contrast to the [0/90]^ 
case, the curves for 0 = 45° and 90° do not switch signs near the hole 
boundary. Also, they become vanishingly small within one ply thickness from 
the hole. 

Interlaminar Shear Stress, a^g--Figs. 9(a) and 9(b) show the radial vari- 
ation of a o for [0/90]^ and [90/0]^ laminates, respectively. On the planes 
of symmetry 0=0° and 90°, p^q is zero. The p^q curves have very steep 
gradients close to the hole boundary, suggesting the presence of a singularity 
there. At two elements away ((r - R)/h = 0.002) from the boundary, the ^Z0 

was about 60 percent as large as S . The two distributions in Figs. 9(a) 

y 

and 9(b) are identical except for their signs. This sign difference was 
expected, as discussed in the previous section. 
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Interlaminar Shear Stress, Figs. 10(a) and 10(b) show the radial 
variation of interlaminar distributions. In both figures the 
curves for 6=0° and 90° do not approach zero at the hole boundary, as might 
be expected. On the hole boundary, of course, the corresponding complementary 
stress 0 ^^ should be zero because this boundary is to be stress free. This 
apparent discrepancy can be explained by the presence of an interlaminar stress 
singularity at the hole boundary. As shown in Ref. 7, a shear stress and its 
complement (e.g., and a^^^need not be equal at a stress singularity. 

Circumferential Distributions 

This section presents the circumferential distributions for a„, a , and 

0 z 

o^g near the hole boundary. Again the stresses from the two elements closest 
to the hole were neglected. Therefore, the stresses shown closest to the hole 
are at r/R = 1.0001, that is, at a distance of 0.002h from the hole boundary. 

Hoop Stress, ag--Fig. 11 shows the Og distribution around the hole for 
both laminates. In this figure ply-thickness average values of Og are plotted 
for each ply. The ply stresses calculated from the 2-D exact solution are 
included for comparison. The ply-thickness average Og in the 0° ply agrees 
very well with the exact solution in the region 0° < 6 < 45°, Similar agree- 
ment is shown for the 90° ply in the region 45° < 0 < 90°. The largest dif- 
ference between the two solutions was about 15 percent. 

Interlaminar Normal Stress, a^--Fig. 12 presents the interlaminar 
distribution around the hole. The symbols in this figure are stress esti- 
mates that will be discussed in the next section. For the [0/90]^ laminate 
(Fig. 12(a)), is compressive for most of the region along the hole 
boundary. The stress is tensile only in the regions 10° < 0 < 28° and 
80° < 0 < 90°. Its largest compressive value occurs around 0 = 60°. 
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For the [90/0]^ laminate (Fig. 12(b)), is again compressive for most 
of the region around the hole, with a small tensile region near 0 =90°. The 
largest compressive value occurs near 6 = 60°. 

Interlaminar Shear Stress, interlaminar distributions 

around the hole are shown in Figs. 13(a) and 13(b). Again, the symbols are 
estimates to be discussed in the next section. Except for different signs, the 
distributions are identical for the two laminates (compare Figs. 13(a) 
and 13(b)). These distributions both have their largest values at about 75° 
from the load axis. This suggests that the a^g stress singularity is proba- 
bly strongest near -6 = 75°. Also, the largest a^g value in Fig. 13 is about 
1.65 Sg, which is about 10 times as large as the largest or value 

computed for the same distance (r/R = 1.0001) from the hole.- This comparison 
suggests that the a^g stress singularity is stronger than those for the 
and 0 ^^ interlaminar stresses. 

Quasi -3-D Estimates 

As mentioned earlier, the interlaminar a, and a,, stresses were also 

z zu 

estimated using the quasi-3-D procedure described in the Appendix. This pro- 
cedure was based on the assumption that the curved hole boundary could be 
approximated by a series of straight-edge specimens. Each such segment was 
analyzed by a quasi -3-D method from the literature (e.g., see Ref. 3). The 
tangential strain Cg in each segment was calculated from the exact 2-D 
solution for a laminate hole, as presented in Ref. 1. As explained in the 
Appendix, this Eg strain was applied to the straight-edge segment in two 
different ways. For the simpler case, the Eg value at the hole boundary was 
applied uniformly across the width (radial direction) of each segment. For the 
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more refined second case, each segment was subjected to the nonuniform radial 
distribution of eg obtained from the exact' 2-D solution in Ref. 1. 

Interlaminar Normal Stress, a^— The symbols in Fig. 12 show the 
estimates at 15° intervals along the hole boundary. The square symbols repre- 
sent the results from the uniform eg approach and the circular symbols 
represent the nonuniform Sg approach. For the [0/90]^ laminate (Fig. 12(a)), 
these two approaches estimate the values very well except at 6 = 90°, 
with both approaches giving very nearly the same values. This leads to the 
conclusion that the radial gradient of the £g distribution (see the Appendix 
and Fig. 15) had very little influence on the results. For the [90/0]^ lami- 
nate (Fig. 12(b)), the estimates are also very good. The small discrepancies 
between the estimates and the 3-D finite-element values may be caused by the 
strain gradients in the 6 -direction which were not accounted for by the 
straight-edge analyses. 

Interlaminar Shear Stress, <^ 26 “"^^^’ presents the a^g estimates. 

For both laminates the predictions are in excellent agreement with the 3-D 
values; only the peak value at 0 = 75° was overestimated. Again, both the 
uniform and nonuniform strain approaches yielded nearly identical values. 

This suggests that the Cg gradients away from the hole in the nonuniform 
approach had very little influence on the estimates and that the simpler uni- 
form Cg approach can be used to estimate and o^g around the hole. 

Through-the-thickness and radial distributions of and a^g were 
also estimated using the quasi-3-D procedure. Unfortunately, the agreement 
between these estimates and the corresponding 3-D results was not as good as 
that shown in Figs. 12 and 13. The quasi-3-D procedure needs further evalua- 
tion, but this is beyond the scope of the present study. 
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CONCLUDING REMARKS 

Stress distributions were calculated in [0/90] and [90/0] laminates 
with a central circular hole, using a three-dimensional finite-element analy- 
sis. The laminate was idealized using 20-node isoparametric brick elements. 
The stress distributions were presented in three ways: through the thickness 
at the hole boundary, along radial lines at the interface, and around the hole 
at the interface. 

The interlaminar normal stress a, and shear stress a,, between the 

Z Zv 

0/90 and 90/0 plies showed very steep gradients very near the hole boundary. 
These results suggest that an interlaminar stress singularity exists at the ; 
hole boundary. For both [0/90]^ and [90/0]^ laminates, the calculated inter- 
laminar stresses were similar, with only minor differences. 

The singularity appears to be strongest at about 60° from the load 
axis. But this interlaminar was compressive over most of the region 
around the hole. The o^q singularity appears to be strongest at about 75° 
from the load axis, and it appears to be stronger than either the or 
singularities. The o^g stresses calculated very near the hole boundary were 
about 10 times the corresponding and values. 

A simple procedure was proposed to estimate and a^g at the inter- 
face very near the hole boundary. This procedure was based on a quasi-three- 
dimensional straight-edge model, subjected to the tangential hole-boundary 
strain computed from an exact 2-D analysis. The estimates agreed very well 
with the computed 3-D stresses. Furthermore, these and a^g inter- 
laminar stresses appear to be governed by the tangential strain value at the 
hole boundary, rather than its radial distribution. 
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The simple procedure introduced in this paper appears to be an attractive 
and economical method for estimating the interlaminar stresses very near holes 
in composite laminates. This simple procedure may prove to be a viable alter- 
native to a complete three-dimensional analysis of a laminated plate contain- 
ing a hole. However, the range of applicability of this procedure needs to be 
checked for more complex laminates than those studied in this paper. 
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APPENDIX 

STRAIGHT-EDGE APPROXIMATION FOR A LAMINATE WITH A, HOLE 

The purpose of this Appendix is to describe a simplified quasi-3-D proce- 
dure to obtain and values very near the hole boundary (r/R = 1.0001 

in the present analysis). 

Fig. 14 shows one-eighth of the laminate from Fig. 1. The hole boundary 
was approximated by seven discrete edge segments as shown. Each of these seven 
segments was analyzed as though it were a laminate with a straight free edge, 
using the quasi-3-D procedure described in Refs. 3, 5, and 6. According to 
this procedure, each edge segment was analyzed using a finite-element model for 
the plane normal to the edge. As shown schematically in Fig. 14, such straight- 
edge finite-element models were in radial planes at 15° intervals around the 
hole. These straight-edge models, which are like the model shown earlier in 
Fig. 3, are described in detail in Ref. 3. 

For the [0/90]^ laminates with a hole, the following stacking sequences 
were used in the straight-edge models: [90/0]^, [75/-15]^, [60/-30]g, [±45]^, 

[30/-60]g, [15/-75]g, and [0/90]^. These seven stacking sequences may be 
expressed concisely as [(90 - e)/- 6]5, where 6 is measured from the x-axis. 

For the [90/0]^ laminate with a hole, the corresponding expression is 
[-6/(90 - 0)]g* These stacking sequences were also used in the straight-edge 
analyses in Ref. 3. 

As shown schematically in Fig. 14, a quasi-3-D, straight-edge model is 
subjected to a strain normal to the plane of the model. For the present models 
in radial planes, this strain corresponds to the tangential strains Cq around 
the laminate hole. The tangential strains Eq along radial planes were 
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computed using the 2-D exact solution of Ref. 1. These e„ distributions are 

U 

shown in Fig. 15. 

Two different approaches were used for the strain inputs in the quasi-3-D 
analyses. In the first approach, a uniform Eq was imposed. The uniform 
strain equaled the Eq value at the hole boundary. (For example, for the 90° 
radial plane, a uniform Eg of 4.39 x 10"^ was used.) In the second 
approach, nonuniform distributions of Eg from Fig. 15 were imposed, as shown 
schematically in Fig. 14. 

The interlaminar stresses and a^g were computed (at r/R = 1.0001) 
for each radial plane with both uniform and nonuniform Eg distributions as 
inputs. The results from these quasi-3-D analyses are presented in Figs. 12 
and 13. 
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Figure 1.- Laminate configuration, loading, and stresses. 
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Figure 5.- cTqq distributions through the laminate thickness 
at the hole boundary (r/R = 1). 
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Figure 6.- cr^ distributions through the laminate thickness 
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Figure 12.- cr^ distributions around hole (at the ply 

interface^ z = h). 
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Figure 13.- cf^q distributions around hole (at the ply 

Interface^ z = h), 
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Figure 14.- Straight-edge approximation for quasi-3-D analysis of 

laminate hole. 
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